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We provide a self-contained account of the B vs. H constitutive rela- 
tion near Hd in Type II superconductors with various types of quenched 
random disorder. The traditional Abrikosov result B ~ [ln{H — Hd)]^'^, 
valid in the absence of disorder and thermal fluctuations, changes signifi- 
cantly in the presence of disorder. Moreover, the constitutive relations will 
depend strongly on the type of disorder. In the presence of point disorder, 
B ~ {H — Hci)^^"^ in three-dimensional (thick) superconductors, as shown by 
Nattermann and Lipowsky. In two-dimensional (thin film) superconductors 
with point disorder, B ~ [H — Hd)- In the presence of parallel columnar 
disorder, we find that B ~ exp[—C3/{H — Hd)] in three dimensions, while 
B ~ exp[—C2/{H — ffci)^/^] in two dimensions. In the presence of nearly 
isotropically splayed disorder, we find that B {H - Hdf^^ in both two and 
three dimensions. 



I. INTRODUCTION 

The physics of vortex lines in high-temperature superconductors has attracted much ex- 
perimental and theoretical interest The competition between interactions, pinning, and 
thermal fluctuations gives rise to a wide range of novel phenomena that are both interesting 
in their own right and technologically important. Here, we focus primarily on the effects 
of disorder on vortex behavior, phenomena which are also important for low-temperature 
Type II superconductors. In a disorder-free sample, vortex lines always flow in response to 
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a current, leading to a resistance even at arbitrarily small currents. By contrast, defects 
in superconductors attract vortex lines. Just as a few nails can hold a carpet in place, an 
entire vortex line system can be held in place by a few defects, leading to a large critical 
current — provided the vortices are in a solid phase, where they are held in place relative 
to each other by their mutual repulsion, which gives rise to a shear modulus. If the tem- 
perature is raised so that the vortices are in a liquid state, then vortices held in place by 
the disorder are pinned, but the other vortices can flow around them, leading again to zero 
critical current. 

In the case of the high-temperature superconductors, disorder is naturally present in 
the form of oxygen vacancies. The quantity of such vacancies can be altered by changing 
the doping of the crystal, i.e., by introducing more or less oxygen during the growth pro- 
cess. Such changes can have a strong impact on the pinning of the vortices, and hence 
on the critical current. Substitutional defects can play a similar role in low-temperature 
superconductors. 

Artificial defects added to superconductors are particularly effective in increasing the 
critical current. Of special note are columnar defects, which are generated by bombarding 
the sample with heavy ions that produce damage tracks in their wake These tracks 

pin vortices strongly because their width (~ 60A) is comparable to the vortex core size 
(~ 20 A). When the columnar pins are aligned with the direction of the magnetic field (and 
hence with the vortex lines) a large increase in the critical current is observed. 

However, an even greater increase in the critical current is observed when the columns 
are splayed, i.e., they are not all oriented in the same direction The superior pinning 

properties of columnar defects with controlled splay was predicted theoretically based on 
expected reduction in a variable-range hopping vortex transport mechanism and enhanced 
vortex entanglement due to splay p!0| , p!l| . Despite the technological importance, much re- 
mains to be done in order to understand the behavior of vortex lines in the presence of 
splayed columnar disorder. 

Even at high temperatures, when the analysis is expected to be most tractable, stan- 
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dard approaches seem to break down and give nonsensical results. For example, the boson 
mapping ||T2|-p^] works quite well in describing many of the properties of vortex lines in the 
presence of point or unsplayed columnar disorder. This approach utilizes a formal corre- 
spondence between vortex trajectories and the world lines of fictitious quantum-mechanical 
bosons in two dimensions. In this analogy, the temperature T plays the role of Planck's 
constant the bending energy (or line tension) g plays the role of the boson mass m, and 
the length of the sample L corresponds to jSU for the bosons. The analogy works best in the 
L oo limit, which corresponds to the T ^ limit for the bosons. In this limit, in the 
absence of disorder and interactions, the bosons should form a condensate. In the presence 
of disorder or interactions, some "bosons" are kicked out of the condensate, resulting in a 
condensate density uq which is less than the total density n. For vortex lines, the "conden- 



sate density" is a measure of the degree of entanglement ||T2|,|T3|,|T5[ . This condensate density 
can be calculated for the various types of disorder discussed above. It is well-behaved for 
point and unsplayed columnar disorder, but in the presence of even weak splayed columnar 
disorder, diverges |T^, suggesting that the boson mapping is flawed in the presence of 
splayed columnar disorder. Indeed, Tauber and Nelson conclude that the super-diffusive 
wandering of the flux lines causes the mapping onto non-relativistic bosons to break down. 
Unfortunately, this mapping — which has provided many insights for the cases of no disorder, 
point disorder, and unsplayed columnar defects — seems to be less suitable for understanding 
the behavior of vortex lines in the presence of splayed columnar disorder. 

Here, we focus on the behavior of vortex lines near i^ci, where the lines are dilute. In 
particular, we predict the B{H) constitutive relation for vortex lines in the presence of the 
various types of disorder discussed above. We begin by reviewing the traditional Abrikosov 
result, expected to hold in the absence of disorder and thermal fluctuations. Each vortex 
that enters the sample will gain a free energy proportional to {H — Hd) per unit length. 
However, there is an energy cost due to repulsive interactions between any two vortices 
proportional to exp(— r/A) for r ^ A, where r is the distance between the vortices and A 
is the London penetration depth. In the dilute limit, the interactions with nearest neighbors 



will dominate over interactions with more distant neighbors, and the free energy density is 
given by 

/ = -ci{H - H,i)n + C2n^/^e-^'/^, (1.1) 

where n is the density of vortices, and Ci, C2, and C3 are constants that can be determined 
in terms of the vortex parameters |T6[. (We leave them general here to better elucidate the 
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structure of the argument.) Upon minimizing / with respect to n, we obtain 
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The dominant behavior may be obtained by substituting n oc Cg on the right hand side, so 
that the magnetic field (given hj B = rupQ, where 0o is the flux quantum) varies inversely 
as the square of the logarithm of H — Hd- Plugging in the relevant parameters ci, C2, and 



C3 for a triangular lattice, one finds 



n 
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1.3) 



_47rA2 {H-Hci)_ 

In the presence of disorder, Eq. (|1.3|) will be modified. With point disorder, we obtain 



B (H — Hci) (as calculated |T8|,|T9[ and measured by BoUe et al. |[T^) in 1+1 dimensional 
samples, where the vortices only have one direction transverse to the magnetic field in 
which they can wander (see Fig. |l]). In 2+1 dimensional samples (Fig. |]), we obtain B ~ 
(H — HciY^'^, in agreement with calculations done by Nattermann and Lipowsky [^. In the 



presence of splayed columnar defects, we find that B ~ [H—HdY^^ both for 1+1 dimensional 
(Fig. 1^) and 2+1 dimensional (Fig. ^ samples. For columnar disorder with unbound disorder 
strength, Larkin and Vinokur argue that B ~ ^c(h-Hci) 2+I dimensions. However, 
we show that for the more physical case of bounded disorder, B ~ exp[—C3/{H — Hd)] in 
2+1 dimensions (Fig. |^), while B ~ exp [-C2/{H - HdY^^] in 1+1 dimensions (Fig. |). In 
addition to these relations, we also estimate the prefactors (up to factors of order unity) 
in terms of physical constants such as the temperature, the disorder strength, and the 
superconducting coherence length. These prefactors are important for comparisons with 
experiment. 
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FIG. 1. A vortex line in a 1+1 dimensional sample of width W in the presence of point disorder. 




FIG. 2. A vortex line in a 2+1 dimensional sample in the presence of point disorder. 




FIG. 3. A vortex line in a 1+1 dimensional sample of width W in the presence of splayed 
columnar disorder. 
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FIG. 4. A vortex line in a 2+1 dimensional sample in the presence of splayed columnar disorder. 
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FIG. 5. A vortex line in a 1+1 dimensional sample of width W in the presence of parallel 
columnar disorder. 




FIG. 6. A vortex line in a 2+1 dimensional sample in the presence of parallel columnar disorder. 



A recent torsional oscillator experiment carried out by BoUe tt al. confirms that 
B ~ {H — Hci) for vortex lines in two dimensions in the presence of point disorder. The 
experiment was done by attaching a thin sheet of superconducting 2H-NbSe2 to a high-Q 
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micro-electromechanical device, and measuring the resonance frequency of torsional oscilla- 
tions of the sample. Because a sample with a magnetic moment M in a field H exerts an 
additional torque r = M x H on the oscillator, the magnetic field which penetrates and 
becomes trapped in the sample can be probed by finding shifts in the resonance frequency 
as a function of applied field H. Small jumps in the frequency as a function of applied field 
were observed which were attributed to individual vortices entering the sample. By 



counting the number of jumps as a function of these authors obtained B ~ {H — Hd). In 
the presence of parallel or splayed columnar disorder, the B vs. H curve could be measured 
by a very similar experiment, where the sample was first irradiated isotropically by heavy 
ions to produce the required disorder. In three dimensions, a similar experiment could be 
done using a long thin needle-shaped sample to avoid demagnetizing effects. For point dis- 
order, this experiment was done in the 1960's with ambiguous results In particular. 



work by Finnemore et al. on Niobium samples seems to indicate B ~ {H — Hd)^, with 
X > 1. However, the decades-old data is too rough to allow quantitative comparison with 
the prediction for point disorder B ~ (H — Hd)^^"^. 

We approach these problems differently for the different types of disorder. Parallel colum- 
nar disorder localizes vortices into finite transverse regions of the sample in both 1+1 and 
2+1 dimensions. Therefore, at low densities, the vortices can easily avoid paying a large cost 
associated with their mutual repulsion if they are located in different areas of the sample. 
In the language of the boson mapping described above, we can approximate the intervortex 
interactions as merely restricting the occupancy of any given localized state to precisely one 
boson [0,^. The form of the B{H) constitutive relation is then determined by the low 



energy tails of the density of states. This idea is explored further in Sec. [Ill D| . 

By contrast, as discussed in Sec. y, in the presence of point disorder or splayed columnar 
disorder, noninteracting vortices will be in extended states. As such, repeated collisions 
between vortices play a crucial role in determining the B{H) constitutive relation, as the 
intervortex interactions attempt to localize each vortex into a cage surrounded by their 
neighbors. The energy lost due to restricted vortex wandering will determine B{H). We 



investigate this case more fully beginning in Sec. || with the problem of a single flux line 
superimposed on a background of disorder. The partition function describing a flux line can 
be mapped onto the noisy Burger's, or KPZ, equation describing the fluctuations of an elastic 
interface in the presence of a random space and time dependent potential influencing the 
interface's progress pH|. This problem has been studied extensively for random potentials 



uncorrelated in space and time (appropriate to point disorder for the vortex lines) p6| , ^ , ^ , 
but more recently has also been investigated for correlated potentials |]28|- pO|] . Provided we 
restrict ourselves to the case where the splayed columnar defects are nearly isotropically 
oriented, we can use previous results for the wandering exponent (, which describes how far 
the vortex line wanders transverse to the magnetic field in a distance / in the field direction by 
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the formula {[r(/) -r(0)]2} where r(z) labels the transverse coordinates of the flux 

line. Columnar defects with an approximately isotropic distribution of splay can be created 



by using neutrons or protons to trigger fission of, e.g., bismuth nuclei in BSCCO []31 |. 

In Sec. pn] , we develop an approximate theory for a dilute set of vortex lines in the 
presence of point and splayed columnar disorder. We first derive the scaling of B vs. H, and 
then focus on the prefactor of this relation. In some experiments the columns do not 

traverse the entire sample, so we also consider the effect that finite column size will have 
on our predictions. Specifically, we expect a crossover from the point disorder behavior at 
extremely weak H — Hd to the splayed columnar behavior at somewhat larger H — Hd. 

II. SINGLE VORTEX LINE 
A. Model 

We begin with a model free energy for a single vortex line in a d dimensional sample of 
thickness L. We label the direction of the magnetic field by r, and the transverse position 
of the line at r by a d — 1 dimensional vector r(r). The free energy then reads |[l|,p!3[|: 



F[rir)] = \gj^ rfr + ^ V[r(r), r]dr, (2.1) 



where g is the hne tension. The pinning potential V^[r, r] arises due to the interaction with, 
say, point disorder or splayed columnar defects. Its mean value merely affects the average 
field Hci at which vortex lines will penetrate the sample, so we subtract it out, and assume 



that V[r, r] = 0. We further assume that the noise is Gaussian with a correlator 



A(r - r', r - r') = V[r, T]V[r', t']. (2.2) 

We focus on the case of "nearly isotropic" splay. This is the defect correlator for a set 
of randomly tilted columnar pins, each described by a trajectory r(r) = R + vr with a 
Gaussian distribution of the tilts v, P[v] oc e~'"'^^^'"o in the limit vo ^ oo [|10|- For nearly 
isotropic splay, the Fourier transform of the correlator, defined by 

r-L 
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A(k,cu) = J ci'^-^r^^rfrA(r,r)e-*(''-™), (2.3) 



is given by A/A; |[T0| , p!5[| , which differs from the truly isotropic limit A/ (A;^ + ct;^)^/^. However, 
using the correlator for nearly isotropic splay simplifies the calculations significantly. In the 
language of the noisy Burger's (or KPZ) equation, truly isotropic splay has both spatial and 
temporal correlations, while nearly isotropic splay has only spatial correlations. By focusing 
on the nearly isotropic splay, we can therefore rely on work done on the KPZ equation 
with spatially correlated disorder and avoid the more complicated case of temporal disorder. 
Moreover, we believe that this difference should not affect the physical implications at large 
length scales significantly [p2 |. 



The partition function associated with the free energy of Eq. ( p.l|) , 

Z(r,r)=y'l?r'(r')e-^['-'(^')]/^, (2.4) 

i.e., the path integral of the Boltzmann factor e~^^^'^'^'^^^'^ over all vortex trajectories r'(r') 
running from r' = to position r(r) at r' = r, obeys the differential equation 

^5^(r,r) 



gv^ + nr,r) 



Z(r,r). (2.5) 



dr 

Eq. ( |2.5D can be further transformed by means of the Cole-Hopf transformation (similar to 
the WKB transformation in quantum mechanics) 
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Z(r,r) = exp [$(r,r)] (2.6) 

into 

= — + — (V$)^ + y(r, r). (2.7) 
or 2g 2g 



This equation, known as the KPZ or noisy Burger's equation has been studied in 
the case of uncorrelated (point-hke) disorder in a great variety of contexts. It was studied 
by Forster, Nelson, and Stephen as a model for the velocity fluctuations in a randomly 
stirred, {d — l)-dimensional turbulent fluid. Later, it reappeared as a model for surface 
roughening proposed by Kardar, Parisi, and Zhang [^. In this context, it has inspired a 
great deal of work, summarized in a review article by Halpin-Healy and Zhang p5| . 



There have also been investigations of the properties of this equation with correlated 
disorder A(k) = A/A;^^ in the context of surface roughening. However, to date, there has 
not been any reliable characterization of the exponents for the entire d vs. p space. Medina, 



et al. have found results that seem to be accurate in two dimensions |Q. Halpin- 
Healy [^] has obtained results that are accurate in two dimensions, and may describe the 
behavior for sufficiently small p in higher dimensions as well. Recently, Frey, et al. have 
found exact exponent solutions not only in two dimensions but also for sufficiently large 
powers of p in higher dimensions. (It is by reference to these exact results that we judge 
the accuracy of Medina, et al. and Halpin-Healy's results.) Although they are not able to 
obtain exact results for all d and p, Frey, et al. make a plausible conjecture for the behavior 
at smaller values of p that is in agreement with Halpin-Healy's short range results and with 
numerical simulations ||3^ . 



Our purpose here is to adapt this work to vortex lines in the presence of splayed columnar 
disorder, and explain some of the consequences for experiments. The application of the 
results for point disorder to experiments is also reviewed. In Sec. [11 B| , we set the stage 
for a more detailed discussion with a simple, self-contained renormalization group treatment 
that gives results in agreement with Frey, et al.'s exact results in its range of applicability. 



In Sec. |II Q we discuss the limitations of this approach. 
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B. Renormalization group 



1. Scaling 

To understand the vortex wandering at long length scales, we use the renormalization 
group method of Ref. |^ to integrate out the short- distance behavior. Since, under renor- 
malization, the coefficients of the term and the (V$)^ will no longer remain identical, 
we replace Eq. ( |2.7D by the more conventional notation 

^ = ^V^$ + ^(V$)2 + I/(r,r), (2.8) 

where initially v = 2\ = T /2g, and we have absorbed a factor of T into V so that A(k) = 
A/{T'^k). We generalize to correlated disorder A(k) = A/{T'^k'^P). It is instructive to 
generalize the correlator this way, because then not only does p = 1/2 describe nearly 
isotropic splayed columnar disorder, but p = generates the results for point disorder as 
well. This renormalization group will be based around a perturbation series in powers of A. 
We ffist rescale Eq. ( p. 81) by a scale factor b, with 

r ^ br, (2.9) 
r ¥t, (2.10) 

(2.11) 

where z and x eventually be chosen to keep various coupling constants fixed under the 
renormalization procedure. Since z describes the ratio of the scaling in the timelike direction 
with that in the spacelike direction, ( = 1/z is exactly the wandering exponent that we are 
looking for. Upon inserting these transformations into Eq. ( p.8| ), we see that the equation 
remains invariant under these changes provided we rescale z/. A, and A via 



u ¥-^u, (2.12) 

A ^ b^+'-^X, (2.13) 

A ^ 52p+i-d+^-2xA. (2.14) 
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In the absence of the nonhnearity {i.e., X = 0), the equation becomes completely scale 
invariant if we choose z = 2 and x = (2p + 3 — d)/2. However, any small A then rescales 
according to 

^ ^ ^(2p+3-d)/2^^ (2.15) 

and thus the nonlinear term will be relevant for d < 2p + 3 {i.e., d < 4 for splayed columnar 
disorder and c/ < 3 for point disorder). For > 2p + 3, we expect the mean field exponents 
displayed in Eq. ( 2.12 - 2!T^ to be accurate, and ( = 1/2. However, in the physical situation 



of interest, d < 2p + 3 and we expect the scaling exponents to change due to the nonhnearity. 

2. Perturbation theory 

To understand the case d < 2p + 3, we rewrite Eq. ( |2.8| ) in Fourier space via the definition 
^k,tu) = J d'^-h £ dT^r,T)e-'^^-''-'^^\ (2.16) 

obtaining 

$(k, iu) = Go{k, iu)V{k, Lj) - -Go{k, lu) [ d'^-^k' dcu'k' ■ (k - k')$(k', cj')$(k -k',iu- lu') 

2 J Jo 

(2.17) 

where 

Go{k,uj)= (2.18) 

We define a renormalized Green's function G{k, u) via 

$(k, uj) = G{k, uj)V{k, uj) (2.19) 

and calculate G{k,u!) perturbatively in A. The perturbation series can be summarized dia- 
grammatically by giving graphical representations to the renormalized and bare propagators 
G{k, u) and Go{k, u), the disorder V{k, u), the interaction, and the disorder correlator A(k) 
as shown in Fig. 0. As expressed in Eq. ( |2.19|) , ^{k,uj) is represented by a double arrow 
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followed by a cross to represent the disorder. Eq. ( p.lTD can then be represented diagram- 
matically as in Fig. We obtain an iterative solution by substituting Eq. ( p. 17] ) for each 
of the $(k, cu) terms appearing on the right hand side of Eq. ( |2.17| ), thereby obtaining an 
iterative series in powers of A. The result for the renormalized propagator is presented 
diagrammatically (to second order in A) in Fig. |^. 



13 



(a) (b) (c) (d) ^ (e) 

FIG. 7. The definitions of the graphical elements of our perturbation series, (a) represents 
the renormalized propagator G(k,uj), (b) represents the bare propagator Go(k, w), (c) represents 
the noise V(k,uj), (d) represents the interaction — ^/q • (k — q), and (e) represents the disorder 
correlator 2A(k). 



+ 




FIG. 8. Diagrammatic representation of Eq. ( 2.17] ) 






FIG. 9. Iterated version of Eq. (2.17) to second order in A 
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We wish to use this perturbation theory to calculate renormalized versions of the pa- 
rameters z/, A, and A (which we will denote by z/, A, and A respectively). We define v 
by 

1 



limG(k,0) ^ 



(2.20) 



Upon multiplying Eq. ( p. IT] ) (or its diagrammatic equivalent, Fig. ||) by V^(— k, —iS) and 
averaging over the noise we obtain the equation for the renormalized propagator rep- 



resented diagrammatically to one loop order in Fig. |T0 




+ 4 



FIG. 10. Diagrammatic expansion of the renormalized propagator, to the one-loop order. 



We define A via a renormalized vertex that contains the effects of the interactions, shown 
on the left hand side of Fig. |TT]. The vertex amplitude, in the limit of small k, q, cj, and 
is given by 

A 



-q • (k - q)Go(k, a;)Go(q, f^)Go(k - q, u; - fi), 



(2.21) 



and serves as the definition of A. Expanding in terms of the bare quantities, we obtain the 



equation determining A to one loop, shown in Fig. 11. 
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FIG. 11. Diagrammatic expansion of the renormalized vertex, to the one-loop order. 



We define the renormalized noise correlator A by 



^*{k,Lu)^{k,Lu) = 2A{k)G{k,uj)G{-k,-Lu) 



(2.22) 



in the limit of k — and u; — > 0, where 



A(k) 



(2.23) 



Expanding in terms of the bare quantities gives rise to Fig. 12 



+ 4 




FIG. 12. Diagrammatic expansion of the renormalized noise correlator, to the one-loop order. 
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For the case p = 0, diagrams like these are expressed in integral form and evaluated in 
detail in Ref. and, e.g., by Barabasi and Stanley |3^. The case p 7^ does not produce 
any new complications, so we simply report the results: 



V = V 

A = A 
A = A 



\ A^A - 2p - 3 /-A ^_ ■ 



1 + 



d-4 



(2.24) 
(2.25) 
(2.26) 



where A is a cutoff in momentum space and is the surface area of a d-dimensional sphere 
divided by (27?)"^. The nonlinear coupling A is unrenormalized, as required by Galilean 
invariance p5| , p6| . Moreover, the noise correlator is also unrenormalized for any p > 0; 
the diagram correcting the vertex produces only white noise (point disorder). This can be 
seen by noting that the one-loop diagram that renormalizes the noise correlator (shown in 
Fig. |T2|) is regular as k — > because the momenta passing through the disorder correlator 



remain finite as k — > 0. Since the other diagrams in Fig. ^diverge as only white noise, 
rather than correlated noise, is produced. We will ignore the effect of this white noise for 
now, since the correlated noise is more singular, and return to discuss its effects in Sec. [11 (J 



3. Renormalization group recursion relations 



All corrections from first order perturbation theory in Eqs. ( |2.24|) — ( p.26|) are well be- 
haved for > 2p + 3, as expected from the earlier scaling argument. In lower dimensions, 
the renormalization group procedure resums this divergent perturbation series by integrat- 
ing over modes with high momentum Ae~' < k < A and rescaling the resulting equations 
by k ^ e^'k. Upon combining the scale transformations Eqs. (|2.12|) — ( p.l4|) with the dia- 
grammatic results above we can easily obtain the flow equations: 

A^A d-2p-3 



z-2- 



du 

m 



rv3 A{d-i) 



-K, 



d-i 



(2.27) 
(2.28) 
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dA 



A 



z-2x-d+l + 2p + 5p 



We can express these in a single flow equation for the combination g"^ = ^3^" 



J- 



dl 



3 + 2p-d , 3(3 + 5p,o)(rf-l)-6p-6 
-fif + Kd-ig 



(2.29) 



(2.30) 



2 - ■ 8(rf- 1) 

We look for fixed points of Eq. ( p.30|) . We first review the case p = 0, i.e., point disorder. 
In d = 2, the equation reads 



dg 1 1 3 

—rr = -Q -f^i Q ■, 

dl T 4 ' 



(2.31) 



which has an unstable fixed point (corresponding to no disorder) at (7 = and a stable fixed 
point at (7 = Upon inserting this fixed point value into Eqs. ( p.27|) — (|2.29| ), we see that 

the long- wavelength physics is characterized by z = 3/2 and x = V^? giving a wandering 
exponent of C = = 2/3, in agreement with other results |]26| , |25| , |28| - pO|] . In = 3, by 
contrast, Eq. ( |2.3CI| ) reads 

The fixed point at (7 = is still unstable, but any small disorder flows off to (7 = 00, or 
strong coupling, where our renormalization group is no longer accurate. Thus determination 
of the wandering exponent for point disorder in three dimensions is beyond the scope of this 



method [25,26 



We now turn to the case p > 0. The flow equation for g now reads 

dg 2p + 3-d 3d - {2p + 3) 3 
-dl = 2 ' + 8 id -I) ■ 



(2.33) 



For < 2p + 3, the coefficient of the linear term is positive, while the coefficient of the cubic 
term is negative, leading to an unstable fixed point at (? = and a stable one &i g = 



4(rf-l) 



Eqs. (12:271)— (^1291) now lead to z = {3 + d-2p) /3, x = (2p + 3 - rf) /3, yielding a wandering 
exponent of 

3 



c 



(2.34) 



3 + d- 2p 

For splayed columnar disorder, p = 1/2, and so the wandering exponent is C = 3/4 in two 



dimensions and C = 3/5 in three dimensions [Pq-|30[|. 



C. Discussion 



One issue that has not been addressed is that of white noise. Recall from Sec. LIB 2| that 



the diagram that renormalizes the disorder correlator A(k) does not produce any correlated 
disorder; however, it does produce white noise. Thus, even if point disorder were not initially 
present, it would be produced by the renormalization group. Naively, one would expect that 
correlated disorder would dominate over white noise for any p > since correlated disorder 
is more singular than white noise at k = 0. This is, in fact the case for p sufficiently large; 
however, for small p the white noise will dominate. Frey et al.. PDI have shown that below 



d = 3 + 2p, the renormalization group allows two fixed points — one long-ranged (dominated 
by correlated disorder), and one short-ranged (without correlated disorder). Moreover, they 
argue that the fixed point with the larger dynamic exponent z will be stable, and that if 
the long-ranged fixed point is the stable one, its dynamic exponent is characterized by the 
exact result 

3 + ^ - 2p 

z\r = 5 , (2.35) 



in agreement with the simpler calculation of Sec. [11 B|. In the case d = 2, we can use the 



known result for point disorder z = 3/2 to show that the short-ranged fixed point is stable 
for p < 1/4, while the long-ranged one is stable for p > 1/4. This establishes that for the 
case of splayed columnar disorder in 2 dimensions, the results are unaffected by the white 
noise. 

The situation is less clear in 3 dimensions because the dynamic exponent of the short- 
ranged fixed point, z^r, is not known. Nevertheless, based on the above considerations, we 
expect that there is a curve pc{d) such that for p > Pc{d), the long-ranged fixed point is 
stable, while for p < Pc{d), the short-ranged fixed point is stable. Frey et al.. conjecture 
that Pc{d) = based on the fact that Pc(2) = 1/4, and Pc(5) = 1. (The latter is known 
from the fact that in d = 5,p = 1, this equation corresponds to the Burger's equation 
with non-conserved noise, previously studied by Forster et al. ||26|) This conjecture leads 
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immediately to a result for Zsr-, namely, z^r = which is in agreement with the results 



of Halpin-Healy p9[ and with numerical simulations in c? = 3. In three dimensions, 
this would therefore imply that splayed columnar disorder is at the boundary between the 
regions of stability between the short- and long-ranged fixed points, and hence that ( = 3/5 
for both point and splayed columnar disorder. 



III. DILUTE VORTEX LINES 

In this section, we show how the B vs. H constitutive relation, which can be measured 
experimentally follows from a knowledge of the exponent (. Our central result, that 

Br^{H-H,i)^, (3.1) 

applies whenever the lines are dilute and < C < 1 [^]. However, the prefactor (which 
is important for comparison with experiment) will depend somewhat on the experimental 
regime. Here, we first derive the above scaling relation, generalizing results for point dis- 



order p!8| , |20| , and then find expressions for the prefactor in various regimes of temperature 
and disorder. 

Recall that parallel columnar defects localize the vortices, yielding ( = 0. Because the 
vortices are localized, Eq. ( ^.Ij ) does not apply to this case. In Sec. p.llL)| , we analyze this 
case with the boson mapping, finding 



B ~ exp 



(3.2) 



A. B vs. H scaling relation with point or splayed columnar disorder 

We first review the scaling properties of the free energy. If we rescale the system by 
scaling the r direction by a factor /, then the transverse directions rescale by a factor Z^. 
According to Eq. ( [^.ID , the elastic term of the free energy then scales by a factor 
Because the physics at low temperatures reflects a balance between the pinning and elastic 
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energies, we expect that the pinning term of the free energy scales the same way. Thus, the 
pinning energy on a scale / is given by Up{l) ~ P'^^^. 

The wandering exponent Ar{l) ~ l'^ (where we define Ar(/) = |[r(z + /) — r(2;)]2| ^ ) 
describes the transverse wandering of lines at long scales. At sufficiently high temperatures, 
thermal wandering describes the physics at shorter length scales. We can then match the 
small-scale results onto the large-scale results at the length scales at which both should be 
valid. The exact short length scale mechanism will depend on the experimental conditions, 
so we keep it general for now. We assume that there is a distance in the r direction Ic, a 
transverse distance Xc, and an energy scale Uc above which the results from Sec. || are valid. 
(We will provide expressions for these parameters in Sec |IIIB| .) Then, for / > Ic, 



Axil)=xJp\ , (3.3) 



Upil) = (r) , (3.4) 



which have the correct long-distance behavior and match with the required values at Ic- 

When a finite concentration of vortices enter the sample, their wandering becomes limited 
by intervortex collisions at a length scale given by Ax{l*) = oq, where is the average 
spacing between the vortex lines; i.e., at /* = Idao/xc)^^'^ . To find the optimal density of 
vortex lines, we can balance the energy gain (per unit length) g{H — Hd)/ Hd of allowing 
a vortex line to penetrate with the pinning energy lost (per unit length) Up{l*)/l* due to 



collisions [^. This yields a vortex spacing 

'gUH - Hci) 



ao = Xc 



2{C-1) 

(3.5) 



UcHci 

In this, and subsequent formulae, we neglect dimensionless constants of order unity. The field 
B which penetrates the superconducting sample is related to Oq via B = 0o/('2o~^^^~'^)) 
where for the case of vortex lines confined to a plate-like geometry (as in Figs. |T] and W 
is the width of the sample in the third dimension. It follows that 

00 \ glciH-Hc^) ]^ 
^d-iW^-d [ UcHci \ ■ ^ ' 
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For flux lines in two dimensions, Eq. ( p.6|) reduces to 
H — Hci 



for point disorder {( = 2/3) 
for splayed columnar disorder {C, = 3/4). 



while in three dimensions we have 



Icl 



for either point disorder [pOl or splayed columnar disorder {( = 3/5). 



B. Physics at shorter length scales 

We now estimate the values Ic, Xc, and Uc that appear in Sec. p.11 A|. The physics at 



short length scales, before the effects of disorder build up, is given by application of the 
naive scaling analysis of Sec. [11 B 1| to Eq. ( [^.8| ), which gives 2; = 2, or = 1/2. The physics 
is similar to that of a random walk, dominated by thermal disorder, as a function of the 
time-like paramter /: = ul for I < Ic, where u = T/2g. We expect the system to cross over 
to the large-scale behavior when the corrections to this diffusion term become comparable 
to the initial value. Upon generalizing Eq. ( p.24|) so that we only integrate out to a length 



scale Xc, we see that the criterion which determines Xc is simply 



J^_^ dqq'^-^P-^ ^ 1, (3.9) 



rv3 

where we have neglected factors of order unity. This leads to expressions for the crossover 
parameters 



3+2p-d 

Xc={^\ (3.10) 



3+2p-d 

^c = ^{^] (3.11) 



Uc = T (3.12) 

provided d < 3 + 2p. The last equality results from noting that U is unrenormalized when 
( = 1/2. The inequality d < 3 + 2p is satisfled for splayed columnar disorder in two or three 
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dimensions, and for point disorder in two dimensions, but not for point disorder in three 



dimensions. In the latter case, one finds 



Xc = ^e^, (3.13) 
/e = ^e— , (3.14) 
Uc = T. (3.15) 

The above results apply at sufficiently high temperatures. However, as is evident from 
Eq. (|3.10|) , Xc decreases with decreasing temperature. If Xc < C,, where ^ is the (transverse) 
cutoff provided by the superconducting coherence length, then the thermal regime is absent 
entirely. From Eq. ( p.lO| ), we see that this breakdown occurs for temperatures T < T*, 
where 

T* = (gA^^^"-'')'^'. (3.16) 

Below this temperature, we must use a zero-temperature treatment to determine the 
characteristic scales Xc, Ic, and Uc- Consider the free energy contributions displayed in 
Eq. ( p.l|) , given that the vortex line has typically wandered a transverse distance Xc = ^ 
in a longitudinal distance Ic- We assume for simplicity that ^ <^ Ic- The energy cost 
of this wandering arising from the first term of Eq. {2A) is approximately g^'^/lc- This 



energy is offset by the line's ability to find a more hospitable pinning environment. Let 
V{lc) = V[r{T),T]dT describe the pinning energy of the wandering line. The gain in 
energy due to wandering should be of order the standard deviation of this zero mean random 



variable, namely yV'^{lc)- We have 

VW) = t t dT'V[T{r),r]V[v{T'),T'] (3.17) 
Jo Jo 

= /^A / 7^— r^T^- (3-18) 

The final integral has an infrared cutoff given by ag ^ and an ultraviolet cutoff given by 
due to the finite size of the vortex core, the vortex line only sees a different disorder profile 
when it wanders a distance ^. For d > 1 + 2p, as is the case for point disorder in two 
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or three dimensions and for splayed columnar disorder in three dimensions, the ultraviolet 
cutoff dominates and we have 

LA 



y'i^c) - (3-19) 
For splayed columnar disorder in two dimensions, we find a logarithmic correction. 



V^{h) = lcA\n{ao/0- (3-20) 

Balancing the energy gain due to disorder with the energy loss due to wandering leads 
to a characteristic length Ic 

/ 2f<i+3-2p\ 1/3 

Ic = j (3.21) 

for (i > 1 + 2p, and 

for d = 1 + 2p (splayed columnar disorder in two dimensions). The corresponding energy 
scale is 

except for splayed columnar disorder in 2 dimensions, where 

U,= (gAeH(^o/OY^' (3.24) 



Up to logarithmic corrections, these results match smoothly onto the high-temperature 



formulae of Eqs. (|3.1CI|) — ( p.l2|) in the region where they both apply, namely, T 



The above results require Ic ^ ^, an assumption that breaks down if A is sufficiently 
large. In fact, the results observed for point disorder in two dimensions by Bolle et al. [|1^] in- 
dicate that for this experiment, Ic < ^. The relevant estimates in this regime are summarized 
in Appendix 0. 
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Combining the results for the matching parameters from Eqs. (p.lO|) — ( p.l2|) and ( |3.21|) - 
( p.24| ) with the B vs. H constitutive relation of Eq. ( |3.6| ) leads to 

' 005'^ H - Hci 



B = I 



WA H, 



Hci 



for T > {gAif^ 
for T < {gAif'^ 



for point disorder in two dimensions (Fig. 



B 



¥ K2± fH -H^i\^'^ 

g gA ' 



H, 



cl 



A 



Hci 



for T > (^A)i/3 
for T < (c/A)i/3 



for point disorder in three dimensions (Fig. |^) 

^2 fH-Hci\"^ 



B 



WA V H, 



cl 



WA In 



[h-HcJ 



3/2 



for T > (^A^ 
for T < (^A^ 



2N1/3 



2\l/3 



for splayed columnar disorder in two dimensions (Fig. and 

^2 ^H-Hci\"^ 



B 



A 



H 



cl 



for splayed columnar disorder in three dimensions (Fig. ^). 



(3.25) 



(3.26) 



(3.27) 



(3.28) 



C. Crossover between splayed columnar and point disorder: finite length columns 

Splayed columnar disorder arising from fission fragments often consists of columns with 
a typical length l^oi that is much smaller than the sample size L (as seems to be the case 
in Refs. f^^)- We then expect to observe a crossover from the behavior typical of splayed 
columnar disorder to that of point disorder sufficiently close to Hd- On scales / such that 
^coi I L, the vortex lines feel the finite size of the columns, and thus the behavior 
should be closer to that described by point disorder. However, for /c ^ ^ ^ ^coi, the vortex 
lines behave as if the columns were infinitely long, and thus the behavior is that of splayed 
columnar disorder. In other words, in two dimensions, we expect the B vs. H constitutive 
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relation to be i? ~ {H — Hd) at very weak fields, where the length scale between collisions is 
above /cob and B ~ [H — HdY^^ at somewhat stronger fields. In three dimensions, since the 
constitutive relation is the same for both point and splayed columnar disorder, the crossover 
will appear only in the amplitude of the power law. 

In the regime / 3> /coi, we expect that the behavior can be described via the methods 



of Sec. [Ill A|, with splayed columnar disorder playing the role of the small-scale mechanism 



alluded to near the beginning of Sec. [Ill A| . Specifically, let Xx and Ux be the transverse 



length scale and energy at which the behavior will cross over from splayed columnar to 
point disorder. (These will play the role of Xc and Uc respectively, while lco\ will play the 
role of Ic-) Then, applying Eq. ( [3.3[ ), we find 



Xc(-j , /c < ^ < ^col 



•^X I ~i I 5 ^ ^COI- 



Ax(/) 

^ ^col / 

Matching these formulae at /cob we see that 

3/4 

Similarly, by using Eq. (p.4[) and matching at /cob we obtain 



(3.29) 



Xx = Xc\-^\ . (3.30) 



Ux = t/c (^y^j . (3.31) 



\ 1/2 



Eq. ( ^.7[ ) then leads to 



! ^^^ '^"^ — for sufficiently weak fields 
XxWUx Hci ^ 
^o9^fH-H^\'/\ ^ (3.32) 
57Tr — — tor stronger nelds. 
x^WUp V H,i J 

We need to find what the field strength will be at crossover. Crossover occurs when the 
distance between vortex lines Oq becomes comparable to Xx- In other words, we cross over 
to the splayed columnar disorder result at the field at which a vortex line typically collides 
with another vortex line every /coi in the z direction. This yields 

B = A^ = A^{}A"'\ (3.33) 

XxW X,W\lcol) ^ ' 
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To find the H at which this occurs, we use Eq. (|3.32| ). Whichever expression we use, the 
same result is obtained, which demonstrates the self-consistency of our result, 

H — Hci Uc 



cl 



(3.34) 



col 



In summary, we conclude that 



B = I 



1/4 



H — Hci H — Hci ^ Uc 



XcWUc \lcoi J Hci ' Hci gy/lJcol 
^og'^'ir fH-Hc,\"' H- He, ^ Uc 



(3.35) 



- XcWUc^^ \ Hci J ' Hci gVU^i 
The parameters Xc-, Ic, and Uc appearing in these equations are those of two dimensional 
splayed columnar disorder, which dominates at short length scales. This yields for high 
temperatures (T ^ [gA^'^]^^'^) 



B = < 

. WA V Hci 
and for low temperatures (T <^ [gA^'^Y^^) 



0^3/2 H_He^ H - Hci A^/^ 

,1/2 



WllitA^/^ Hci ' Hci 



^fH-Hci\^/^ H-Hci 

I cl\ Cl ^ 



Hci 



,1/2 



B 



H-H 



cl 



WllitAV^ 



In I ^ ^ ^col Hc\ 
1 33/2^ H-H,^ 



Hci 



for 



H-H 



In 



cl 



Hci 



^ f gel 



3/2' 



1/2 



while 



B 



WA In 



^ gel ^ 

5g2 \H-H,J 



3/2' 



H-Hci \'^^ 
Hci 



for 



Hci 



In 



cl 



> 



V^-gcl 

7V2 

ghoi 



1/2 



(3.36) 



(3.37) 



(3.38) 



(3.39) 



(3.40) 
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Vortex lines 


Bosons 


9 


m 




h 




(3h 


[H - i/,i)0o/47r 




5/00 


n (boson density) 


Vortex lines in three-dimensional samples 


Two-dimensional bosons 


Vortex lines in two-dimensional samples 


One-dimensional bosons 


Parallel columnar disorder 


Point disorder 



TABLE I. Detailed correspondence of the parameters of the vortex line system with the 
parameters of the boson system. 



D. B{H) constitutive relation with parallel columnar disorder 



The boson mapping p^- p!4| , p9| is particularly useful to study vortex lines in the pres- 
ence of parallel columnar defects. The dimensionality of the fictitious bosons is one lower 
than that of the superconducting sample, i.e., vortices in three-dimensional superconductors 
are described by two-dimensional bosons, while those in two-dimensional superconductors 
correspond to one- dimensional bosons. Parallel columnar disorder plays the role of point 
disorder, while the temperature T plays the role of Planck's constant h, the bending energy 
g plays the role of the boson mass m, and the sample length L plays the role of jSh for the 
bosons. (See Table I for a summary.) Since all eigenstates are localized in one and two 
dimensional quantum mechanics, C = fo^' vortex lines in the presence of disordered paral- 
lel columnar pins in both two and three dimensions. Thus, Eq. ( |3.1| ) does not apply: the 
physics leading to the B{H) constitutive relation is very different. Rather than restricting 
vortex wandering, intervortex repulsion will assign an energy cost to two vortex lines that 
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occupy the same localized region. For B ^ Hd, this energy cost will be prohibitive, and 
we approximate the effects of the interaction as prohibiting multiple occupancy of the same 
state [|13|,^. Thus, the vortex interactions play the role of the Pauli exclusion principle. 



and, in this approximation, the behavior is the same as for spinless non-interacting fermions. 

From Table I, the n(/i) relationship at T = for the fictitious bosons yields the B{H) 
relationship in the thermodynamic limit L — > oo. But the n{^) relationship is simply given 
by 

= r g{E)dE, (3.41) 

J —oo 

where g{E) is the non- interacting density of states per unit energy per unit area. Thus the 
B{H) relation in the dilute limit is determined by the low energy tail of the density of states. 

Larkin and Vinokur |^ determined B{H) in the following fashion: they assumed a 
Gaussian disorder potential, 

{V{r)V{r')) = A,S\r-r'), (3.42) 

from which it follows that at low energies [Q, 

(7(i?)~e-2-9^/^°, (3.43) 

with Eq = Aig/T"^ in the vortex line language. The end result is that 

B e^'^°''^~^''^^/^° , (3.44) 

where A/" is a numerical factor. However, we do not believe this to be an accurate description 
of real flux lines. The disorder is taken to be Gaussian, and as such is not bounded below. 
Therefore, there are states at arbitrarily low energy [E —>■ — oo), as Eq. ( |3.43| ) shows, leading 
to vortex penetration at arbitrarily small fields. Indeed, according to Eq. ( |3.44| ), we would 
expect a small density of vortex lines parallel to the z-axis to penetrate the sample in the 
limit Hz = 0, and even for < 0! This unphysical behavior is an artifact of choosing a 
disorder potential that is not bounded below. To fix this problem, we choose the pinning 
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(3.45) 



potential V{r) from a uniform distribution over the range Eq < V < Ei. While the real 
distribution of the disorder will be bounded, it may not be of this form. Therefore, at the 
end of this section, we discuss how our results would be altered by choosing a different 
(bounded) distribution. 

Clearly, g{E) is bounded from below hy E = Eq, yielding Hd = The form of the 

density of states as E Eq from above is determined by the frequency of large, rare regions 
where the disorder potential is always near the bottom of the band To find g{E), 

we (see, e.g., Ref. |^) estimate the probability p{R, S) of finding a sphere of radius R with 
all energies within 6 of the bottom of the band as 

Jq) (-E^i — Eq^ 

where ds is the dimension of the fictitious bosons {(Ib = d — 1) and Iq is the (microscopic) 
transverse distance over which the disorder potential is correlated, i.e., the radius of the 
columnar defects. 

Because, in the boson mapping, the kinetic energy takes the form (see Table I) — |^V^, 
the low-energy eigenstate produced by a such an anomalous region will be given approxi- 
mately by 

E^Eo + 6+^, (3.46) 

where c is a numerical factor of order unity. Therefore, the probability of finding a state 
between energy E and E + dE using a sphere of radius R is given by 

dp{R,6) 



p{R,E) 



35 
yielding 



5=E-Eo-^ 



(3.47) 



p{R, E) ~ exp 



In 



gR^ 



(3.48) 



Jo / \ El — Eq 

Note that from Eq. ( |3.46| ), since 5 > 0,the lower limit R at which it is possible to create a 
state with energy E is 
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R 



(3.49) 



V 9 E-Eo 

Upon optimizing Eq. ( |3.48D with respect to R, we find that (up to logarithmic corrections) 

the maximum occurs at the lower limit of R given by Eq. (p.49|) . Thus, the form of the 

density of states at low energies is given by 

Jrr2 1 \ {<i-i)/'2-] 



g{E) ~ exp 



c'T' 



1 



git E-Eo 



(3.50) 



The logarithmic corrections alluded to above will change the factor of order unity in the 
exponential, and may introduce pre-exponential terms. We do not, however, calculate these 
effects, as the results are not independent of the details of the distribution from which the 
disorder has been drawn. In particular, if the distribution is bounded below and does not 
vanish too quickly at the lower bound, only the numerical coefficient of order unity c' and 
pre-exponential terms will change; the form of density of states will be the same. However, 
if the distribution falls off faster than any power law at its lower bound {e.g., where the 
probability of obtaining an energy V ~ exp[—K/{V — Eq)]), then the exponent {d — l)/2 in 



Eq. (|3.50|) may change as well. Thus, presuming the disorder potential does not fall off too 
fast near the bottom of the band, 

H^, ^ (^-i)/2- 



B{H) ~ exp 
where is a constant of order unity 



H — Hci 



(3.51) 



IV. OUTLOOK 

We expect that the results described above will be valid when the lines are dilute, i.e., 
when B ^ Hd. Here we discuss the outlook for an understanding of vortex lines in disordered 
superconductors when B ^ Hd, in which case the effects of interactions between the vortex 
lines must be taken into account more carefully. 

At high temperatures, generalizations of "hydrodynamic" models described by Marchetti 
and Nelson |^ and by Nelson and Le Doussal fl^ should describe the lines quite effectively 
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in the presence of point, parallel columnar, and splayed columnar disorder. These models 
predict a liquid-like state at high temperatures. In the presence of point and parallel colum- 
nar disorder, there may be phase transitions to glassy vortex states at low temperatures 
both in 2 and in 3 dimensions p4|-^,p^. In the presence of point disorder in 3 dimensions, 
two types of glassy phases are possible. For sufficiently weak disorder, the vortex lines form 
a "Bragg glass" in which dislocations do not proliferate |^7|-p2|. At stronger disorder, dis- 
locations enter the sample. However, it is not yet clear if there is a sharp phase transition 
separating this "glassy" state with dislocations from a high-temperature flux liquid. 

The case of splayed columnar disorder in 2 dimensions with dense lines has been inves- 



tigated by Devereaux, Scalettar, Zimanyi, and Moon |5^. They conclude that there is a 
transition to a "splay glass" phase at low temperatures. The situation is less clear in three 
dimensions. We expect that, in contrast to the case of point disorder, there will be no Bragg 
glass in the presence of splayed columnar disorder: the pinning produced by columns in 
random directions attracting the vortex lines is likely to have a much stronger entangling 
effect on the vortices than point disorder. Since the dislocation-free Bragg glass observed in 
the presence of point disorder is only marginally stable to dislocations P7|-^, we expect the 
analogous system with splayed columnar disorder to be unstable to dislocations (especially 
screw dislocations, which cause entanglement). 
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APPENDIX A: ZERO-TEMPERATURE KINK REGIME 



In this Appendix, we adapt the results of Sec. |I11B| to the case where the effective 
temperature is low, so a zero-temperature approach is applicable, but the effective disorder 
strength is much stronger than the elastic energy tending to produce straight vortex lines. 
This appears to be the regime in which the experiments of [|TH[ are performed. This regime 
is characterized by the inequality Ic < ^. It follows that the free energy of Eq. ( |2.1| ) no longer 
applies, because it relies on an expansion of the line tension 



In this case, however, we can approximate the wandering on short scales as being by nearly 
transverse kinks of distance ^, with an energy cost of gC, rather than gtl/S,- Balancing this 
against the energy gain of pinning of Eqs. ( p.l9|) and (|3.2CI|) , we obtain 



Ic = ^\ (A.2) 

for point disorder in two or three dimensions, and for splayed columnar disorder in three 
dimensions, while for splayed columnar disorder in two dimensions we obtain 

Alog(ao/0 

In either case, we find 

Uc = g^, (A.4) 
in agreement with the results of Ref . llSl . 
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